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Abstract 

We propose to treat the (j> 4 Euclidean theory constructively in a sim- 
pler way. Our method, based on a new kind of "loop vertex expansion", 
no longer requires the painful intermediate tool of cluster and Mayer ex- 
pansions. 

1 Introduction 

Constructive field theory build functions whose Taylor expansion is perturbative 
field theory [TJ [2] . Any formal power series being asymptotic to infinitely many 
smooth functions, perturbative field theory alone does not give any well defined 
mathematical recipe to compute to arbitrary accuracy any physical number, so 
in a deep sense it is no theory at all. 

In field theory "thermodynamic" or infinite volume quantities are expressed 
by connected functions. One main advantage of perturbative field theory is that 
connected functions are simply the sum of the connected Feynman graphs. But 
the expansion diverges because there are too many such graphs. However to 
know connectedness does not require the full knowledge of a Feynman graph 
(with all its loop structure) but only the (classical) notion of a spanning tree in 
it. This remark is at the core of the developments of constructive field theory, 
such as cluster expansions, summarized in the constructive golden rule: 

"Thou shall not know most of the loops, or thou shall diverge!" 

Some time ago Fermionic constructive theory was quite radically simplified. 
It was realized that it is possible to rearrange perturbation theory order by or- 
der by grouping together pieces of Feynman graphs which share a common tree 
[31 2] . This is made easily with the help of a universal combinatoric so-called 
forest formula [5J [S] which once and for all essentially solves the problem that 
a graph can have many spanning trees. Indeed it splits any amplitude of any 
connected graph in a certain number of pieces and attributes them in a " demo- 
cratic" and "positivity preserving" way between all its spanning trees. Of course 
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the possibility for such a rearrangement to lead to convergent resummation of 
Fermionic perturbation theory ultimately stems from the Pauli principle which 
is responsible for analyticity of that expansion in the coupling constant. 

Using this formalism Fermionic theory can now be manipulated at the con- 
structive level almost as easily as at the " perturbative level to all orders". It 
lead to powerful mathematical physics theorems such as for instance those about 
the behavior of interacting Fcrmions in 2 dimensions [3 [9] , and to more ex- 
plicit constructions [TU] of just renormalizable Fermionic field theories such as 
the Gross- Neveu model in two dimensions first built in [TTJ [T^] . 

But bosonic constructive theory remained awfully difficult. To compute the 
thermodynamic functions, until today one needed to introduce two different 
expansions one of top of the other. The first one, based on a discretization of 
space into a lattice of cubes which breaks the natural rotation invariance of the 
theory, is called a cluster expansion. The result is a dilute lattice gas of clusters 
but with a remaining hardcore interaction. Then a second expansion called 
Mayer expansion removes the hardcore interaction. The same tree formula is 
used twice once for the cluster and once for the Mayer expansion^, the breaking 
of rotation invariance to compute rotation invariant quantities seems ad hoc 
and the generalization of this technique to many renormalization group steps is 
considered so difficult that despite courageous attempts towards a better, more 
explicit formalization P31 [H] , it remains until now confined to a small circle of 
experts. 

The bosonic constructive theory cannot be simply rearranged in a convergent 
series order by order as in the Fermionic case, because all graphs at a given 
order have the same sign. Perturbation theory has zero convergence radius for 
bosons. The oscillation which allows resummation (but only e.g. in the Borel 
sense) of the perturbation theory must take place between infinite families of 
graphs of different orders. To explicitly identify such families and rearrange the 
perturbation theory accordingly seemed until now very difficult. The cluster 
and Mayer expansion perform this task but in a very complicated and indirect 
way. 

In this paper we at last identify such infinite families of graphs. They give 
rise to an explicit convergent expansion for the connected functions of bosonic </> 4 
theory, without any lattice and cluster or Mayer expansion. In fact we stumbled 
upon this new method by trying to adapt former cluster expansions to large 
matrix 4 models in order to extend constructive methods to non-commutative 
field theory (see [16) for a recent review). The matrix version is described 
in a separate publication |17| . Hopefully it should allow a non-perturbative 
construction of the (f>* 4 theory on Moyal space K 4 , whose renormalizable version 
was pioneered by Grosse and Wulkenhaar [15] , 

1 It is possible to combine both expansions into a single one 1131 . but the result cannot be 
considered a true simplification. 
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2 The example of the pressure of 4 



We take as first example the construction of the pressure of <f)\ in a renormal- 
ization group (RG) slice. The goal is e.g. to prove its Borel summability in 
the coupling constant uniformly in the slice index, without using any lattice 
(breaking Euclidean invariance) nor any cluster or Mayer expansion. 
The propagator in a RG slice j is e.g. 

Cj(x,y)= e- am2 e~^ x - v ^' ia a^ 2 da<KM 2l e- cM3 \ x ^ (1) 

where M is a constant defining the size of the RG slices, and K and c from now 
on are generic names for inessential constants, respectively large and small. We 
could also use compact support cutoffs in momentum space to define the RG 
slices. 

Consider a local interaction A J <t> A (x)d A x = ATr</) 4 where the trace means 
spatial integration. For the moment assume the coupling A to be real positive 
and small. We decompose the <j> functional integral according to an intermedi- 
ate field as: 

J ^ Cj (0)e- AT ^ = J cM<7)e^ los(1+iH) (2) 

where dv is the ultralocal measure on a with covariance 8{x — y), and H = 
A 1 / 2 DjaDj is an Hermitian operator, with Dj = C X J 2 . 

The pressure is known to be the Borel sum of all the connected vacuum 
graphs with a particular root vertex fixed at the origin. We want to prove this 
through a new method. 

We define the loop verted^ V = — gTrlog(l + iH). This loop vertex can be 
pictured as in the left hand side of Figure [1] The trace means integration over 
a "root" xq. Cyclic invariance means that this root can be moved everywhere 
over the loop. It is also convenient to also introduce an arrow, by convention 
always turning counterclockwise for a +iH convention, and anti-clockwise for a 
complex conjugate loop vertex V = — iTrlog(l — iH). 

We then expand the exponential as Y^n^nf- To compute the connected 
graphs we give a (fictitious) index v, v = l,...,n to all the a fields of a given 
loop vertex V v . This means that we consider n different copies a v of a with a 
degenerate Gaussian measure dv({a v }) whose covariance is < a v cr v > > v = 5{x — 
y). The functional integral over dv(a) is equal to the functional integral over 
dv({a v }). We apply then the forest formula of [5] to test connexions between 
the loop vertices from 1 to n. (The lines of this forest, which join loop vertices 
correspond to former <j> vertices.) 

The logarithm of the partition function logZ(A) at finite volume A is given 
by this formula restricted to trees (like in the Fermionic case [4]), and spatial 
integration restricted to A. The pressure or infinite volume limit of loS |f| A " > is 
given by the same rooted tree formula but with one particular position fixed at 

2 To avoid any confusion with the former </> 4 vertices we shall not omit the word loop. 
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the origin, for instance the position associated to a particular root line £q. More 
precisely: 

Theorem 2.1. 

JSs.^ = EaE{n[/*.]}*<'.**>k- m 

Gt{v,x 1o ) = Y[ J d 4 xtd±y t J dvr({a v },{w}) 



{n[ { <«-»'^s4s)]}n^ < 4 > 



where 



• eac/i Zme ■£ of i/ie iree joins two different vertices V v m and at point 
xi and ye, which are identified through the function 8(xi — ye) (since the 
covariance of a is ultralocal), 

• the sum is over rooted trees over n vertices, which have therefore n — 1 
lines, with root £ , 

• the normalized Gaussian measure dvx({o~ v }, {w}) over the vector field a v 
has covariance 

< a v ,a v i >= S(x — y)w T (v, v', {w}) 

where w T (v, v', {w}) islifv = v', and the infimum of the wg for £ running 
over the unique path from v to v' in T if v ^ v' . This measure is well- 
defined because the matrix w T is positive. 




Figure 1: Loop vertices and a tree on them 



This is indeed the outcome of the universal tree formula of [6] in this case. 
To check it, we need only to move by cyclicity the local root of each loop nearest 
to the global root in the tree. This global root point is chosen for simplicity in 
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formulas above at a particular root line £ , but in fact it could be fixed anywhere 
in an arbitrarily chosen "root loop" , as shown on the right hand side of Figure 
[1] (with all loops oriented counterclockwise). 

But there is an other representation of the same object. A tree on connecting 
loops such as the one shown in the right hand side of Figure Q] can also be drawn 
as a set of dotted lines dividing in a planar way a single loop as in Figure [5J 
Each dotted line carries a S(xi — yi) function which identifies pairs of points 
on the border of the loop joined by the dotted line, and is equipped with a 
coupling constant, because it corresponds to an old 4 vertex. This second 
picture is obtained by turning around the tree. The pressure corresponds to the 
sum over such planar partitions of a single big loop with an arbitrary root point 
fixed at the origin, The corresponding interpolated measure dv can be described 
also very simply in this picture. There is now a a v field copy for every domain 
v inside the big loop, a w parameter for each dotted line, and the covariance 
of two a v and ov fields is the ordinary S function covariance multiplied by a 
weakening parameter which is the infimum of the w parameters of the dotted 
lines one has to cross to go from v to v' . The counterclockwise orientation of 
the big loop corresponds to the +iH convention. 



Figure 2: The big loop representation 

In this new picture we see indeed many loops... but the golden rule is not 
violated. In this new representation it simply translates into 

"Thou shall see only planar (or genus-bounded) structures..." 

(Recall that genus-bounded graphs are not many and don't make perturba- 
tion theory diverge.) 

Let us prove now that the right hand side of formula ([3|) is convergent as 
series in n. 

Theorem 2.2. The series is absolutely convergent for A small enough, and 
the sum is bounded by KM 3 . 
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Proof We shall use the first representation of Figure [T] Consider a loop vertex 
V v of coordination k v in the tree. Let us compute more explicitly the outcome 
of the k v derivatives Yli=i So-(x ) ac t m g on 

V = — Trlog(l + iH) 

which created this loop vertex. 
Consider the operator 

Cj(a)= . D,-. (5) 

Calling x\ the root position for the loop vertex V v , that is the unique position 
from which a path goes to the root of T, the loop vertex factor V v after action 
of the derivatives is 

[JI^LylK = i(-*VA) fe "^nQ(a,x TW , a ; T(l+1) ) (6) 

i=l * r i=l 

where the sum is over all permutations r of [2, ...,k], completed by r(l) = 
r(fc + 1) = 1. 

To bound the integrals over all positions except the root, we need only a 
very simple lemma: 

Lemma 2.1. There exists K such that for any x and any v 

|[C 3 -(ov,)]*»(a;,aO| <K^M^~ 2k ^ Va v . (7) 

Since iH is anti-hermitian we have ||(1 +iH)~ 1 \\ < 1. It is obvious from |T]) 
that IICjII < KM~ 2j , hence \\Dj\\ < KM~ j . We have 

[Ci{o- v )\ k "(x,x) = / dydzD j (x,y)A(y,z)D j {z,x) =< f,Af > (8) 



for / = Dj(x,.) and A = (1 + iH)- 1 [(7,(1 + iH)- 1 ]^" 1 . The norm of the 
operator A is bounded by K k *~ l M~'* i{ - k "- 1 ') . Since ||/|| 2 < KM 2j , the result 
follows. □ 
To bound the dx£ integrals we start from the leaves and insert the bound ([7]) , 
which also means that the multiplication operator Cj(a v )] * (x, x) (diagonal in x 
space) has a norm bounded by K kv M^ i ~' 2kv " uniformly in a . We then progress 
towards the root. By induction, multiplying norms, adding the \{— i\fX) kv 
factors from (J5]) and taking into account the factorials from the sum over the 
permutations r in ^ gives exactly 

Y[ -(K - i)\\ k */ 2 K k *M ij - 2jk ". (9) 

V 

For a tree on n loop vertices J2 V ^« = ~ ■"•) nence — ^^w) = 4n — 

4(n — 1) = 4 so that collecting all dimensional factors we get a M 4 -? global n 
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independent factor as should be the case for vacuum graphs in the </> 4 theory in 
a single RG slice. 

We can now integrate the previous bound over the complicated measure dvT 
and over the {we} parameters. But since our bound is independent of a v , since 
the measure dv(o~) is normalized, and since each we runs from to 1, this does 
not change the result. 

Finally by Caylcy's theorem the sum over trees costs }j W^Tn • ^he n ' 
cancels with the \/n\ of ((3]) and the l/(k v — 1)! exactly cancel the one in ([9]) . 
It remains a geometric series bounded by |M •' (AiT)™ -1 hence convergent for 
small A, and the sum is bounded by K.M^ . □ 

3 Uniform Borel summability 

Rotating to complex A and Taylor expanding out a fixed number of </> 4 vertices 
proves Borel summability in A uniformly in j. 

Definition A family fj of functions is called Borel summable in A uniformly 
in j if 

• Each fj is analytic in a disk Dr = {A|RcA _1 > \/R\; 

• Each fj admits an asymptotic power series ^2f I o J jj c X k fa 8 Taylor series 
at the origin) hence: 

r-l 

fj(X) »<> A ' • ( 10 ) 

k=Q 

such that the bound 

\RrjW\<A jP r r\\X\ r (11) 

holds uniformly in r and A G Dr, for some constant p > independent of 
j and constants Aj > which may depend on j. 

Then every fj is Borel summable [19], i.e. the power series X^fc a i,fe|r con " 
verges for \t\ < -, it defines a function Bj{t) which has an analytic continuation 

in the j independent strip S p = {t\ dist (£,R + ) < -}. Each such function 
satisfies the bound 

\Bj(t)\ < Bje* for<GlR + (12) 

for some constants Bj > which may depend on j. Finally each fj is represented 
by the following absolutely convergent integral: 

1 Z" 00 

fj(X) = - e~*Bj(t)dt forAGCfl. (13) 

^ Jo 

Theorem 3.1. The series for the pressure is uniformly Borel summable with 
respect to the slice index. 
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Proof It is easy to obtain uniform analyticity for Re A > and |A| small enough, 
a region which obviously contains a disk Dr. Indeed all one has to do is to 
reproduce the previous argument but adding that for H Hcrmitian, the operator 
(1 + ie^H)- 1 is bounded by y/2 for \6\ < tt/4. Indeed if tt/4 < Argz < 3tt/4, 
we have |(1 + iz)~ 1 \ < V2. 

Then the uniform bounds (jlip follow from expanding the product of resol- 
vents in ([6]) up to order r — 2(n — 1) in A by an explicit Taylor formula with 
integral remainder followed by explicit Wick contractions. The sum over the 
contractions leads to the p r r\ factor in (fTTj) . □ 

4 Connected functions and their decay 

To obtain the connected functions with external legs we need to add resolvents 
to the initial loop vertices. A resolvent is an operator Cj(o~ r ,x,y). The con- 
nected functions S c (x±, ...,X2 P ) are obtained from the normalized functions by 
the standard procedure. We have the analog of formula (|3|) for these connected 
functions: 

Theorem 4.1. 

S c ( Xl ,...,x 2p ) = EE^E(Ii [ f dw ' [d 4 x e d 4 yA 

it n=l ' T ^ ££T J0 J ) 

[ dv T ({<T v },{<T r },{w})\T[ [S(x t -ye)-z — — r^] f 

p 

) , (14) 

v r—1 

where 

• the sum over ir runs over the pairings of the 2p external variables into 
pairs (x^i),^^)), r= l,...,p, 

• each line I of the tree joins two different loop vertices or resolvents V v ^ 
and V v i(j\ at point xg and yi, which are identified through the function 
5(x£ — yi) because the covariance of a is ultralocal, 

• the sum is over trees joining the n + p loop vertices and resolvents, which 
have therefore n + p — 1 lines, 

• the measure di'T({o' v } 1 {o>}, {w}) over the {a} fields has covariance < 
o~a,o~ a ' >= S(x — y)w T (a, a', {w}) where w T (a, a', {w}) is 1 if a = a' 
(where a, a' € {v},{r}), and the infimum of the wi for i running over 
the unique path from a to a' in T if a ^ a' . This measure is well-defined 
because the matrix w T is positive. 

Now we want to prove not only convergence of this expansion but also scaled 
tree decay between external arguments: 



Theorem 4.2. The series \1J$ is absolutely convergent for A small enough, its 
sum is uniformly Borel summable in X and we have: 

\S c {z u z 2p )\ < {2p)\K p \X\ p - 1 M 2p3 e -cMid( Zl ,...,z 2p ) (15) 

where d(zi, z 2p ) is the length of the shortest tree which connects all the points 
01 j . . . , Zp. 

The proof of convergence (and of uniform Borel summability) is similar to 
the one for the pressure. 

The tree decay (fT5|) is well known and standard to establish through the 
traditional cluster and Mayer expansion. It is due to the existence of a tree 
of Cj propagators between external points in any connected function. In the 
present expansion, this tree is hidden in the resolvents and loop vertices, so that 
an expansion on these resolvents (and loop vertices) is necessary in one form 
or another to prove (|15p . It does not seem to follow from bounds on operator 
norms only: the integral over the a field has to be bounded more carefully. 

The standard procedure to keep resolvent expansions convergent is a so- 
called large/small field expansion on a. In the region where a is small the resol- 
vent expansion converges. In the large field region there are small probabilistic 
factors coming from the dux measure. This is further sketched in subsection 

EES 

However the large/small field expansion again requires a discretization of 
space into a lattice: a battery of large/small field tests is performed, on the 
average of the field a over each cube of the lattice. We prefer to provide a 
new and different proof of (fT5|) . It relies on a single resolvent step followed by 
integration by parts, to establish a Fredholm inequality on the modulus square 
of the 2p point function. From this Fredholm inequality the desired decay follows 
easily. The rest of this section is devoted to the proof of (fT5|) in the simplest 
case p = 1 . The most general case is sketched in subsection 15.11 

The two point function S c is simply called S(x,y) from now on, and for 
p = 1 (fT5|) reduces to 

\S(x,y)\ < KM 2j e~ cMi \ x - y \. (16) 

We work with n, T and {w} fixed in (fT4")l . We use the resolvent as root for T, 
from which grow q subtrees T±, ...,T 9 . In more pictorial terms, (|14p represents 
a chain of resolvents from x to y separated by insertions of q subtrees. Figure [3] 
is therefore the analog of Figure [1] in this context 

A representation similar to the big loop of Figure [2] pictures the decorated 
resolvent as a half-circle going from x to y, together with a set of planar dotted 
lines for the vertices. The +i convention again corresponds to a particular 
orientation. For reason which should become clear below, we picture the planar 
dotted lines all on the same side of the x-y line, hence inside the half- disk. 

To each such drawing, or graph G, there is an associated Gaussian measure 
dvc which is the one from which the drawing came as a tree. Hence it has a field 

3 A similar figure is a starting point for the 1PI expansion of the self-energy in [7] [9]. 
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Figure 4: The half-circle representation of Figure [3] 



copy associated to each planar region of the picture, a weakening parameter w 
associated to each dotted line, and the covariance between the a fields of different 
regions is given by the infimum over the parameters of the dotted lines that one 
has to cross to join these two regions. 

There is also for each such G an amplitude. Let us write simply J dvQ for 

the normalized integral J^Y[eeG^ We I dvc({o~}, {w}). If the graph has n dotted 
lines hence 2n + 1 resolvents from x to y, its amplitude is 

2n+l 

dv G / [Jp 4 ^] II Cite.Si-i.a*) (17) 

J f £ G i=l 

where the product over t runs over the dotted lines and the product over % runs 
over the resolvents along the half-circle, with xq = x and X2n+i = V- &i is the 
field copy of the region just before point Xi and the 2n positions x±, ..., xi n are 
equal in pairs to the n corresponding x^s according to the pairings of the dotted 
lines. 

We shall prove 

Lemma 4.1. There exists some constant K such that for A small enough 

sup \A G (x,y)\ < (|A|iO n/2 M 2 'e- cM3 >-< (18) 

G,n(G)=n 

From this Lemma (fTBf obviously follows. Indeed the remaining sum over 
Cayley trees costs at most K n n\, which is compensated by the ^ in (jT4j) . In 
the language of planar graphs the planar dotted lines cost only K n . Hence 
the sum over n converges for A small enough because of the |A| n / 2 factor in 
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fTS} . Remark that this factor |A|™/ 2 is not optimal; |A|™ is expected; but it is 
convenient to use half of the coupling constants for auxiliary sums below. 

We apply a Schwarz inequality to \Ac(x,y)\ 2 , relatively to the normalized 
measure dvQ\ 

\A G (x,y)\ 2 < Agug( x >v)i (19) 
Agug( x >v) = I dv G j [ Y[ d 4 x e d 4 Xi\ 

J ^ £<EG 

2n+l 

[ Cj{ai,Xi- 1 ,Xi)Cj{a i ,Xi- 1 ,Xi) (20) 
i=i 

with hopefully straightforward notations. 

The quantity on the right hand side is now pointwisc positive for any a . 
It can be considered as the amplitude A GuG (x,y) associated to a mirror graph 
GUG. Such a mirror graph is represented by a full disk, with x and y diametrally 
opposite, and no dotted line crossing the corresponding diameter. The upper 
half-circle represents the complex conjugate of the lower part. Hence the upper 
half-disk is exactly the mirror of the lower half-disk, with orientation reversed, 
see Figure O 




Figure 5: The mirror graph G U G for the graph G of Figure [4] 

The Gaussian measure associated to such a mirror graph remains that of G, 
hence it has a single weakening w parameter for each dotted line and its mirror 
line, and it has a single copy of a a field for each pair made of a region of the disk 
and its mirror region. Let's call such a pair a "mirror region". The covariance 
between two fields belonging to two mirror regions is again the infimum of the 
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w parameters crossed from one region to the other, but e.g. staying entirely in 
the lower half-disk (or the upper half-disk). 

We shall now perform a single resolvent expansion step and integration by 
parts, together with a bound which reproduces an amplitude similar to A Gkj q. 
The problem is that the category of mirror graphs is not exactly stable in this 
operation; this bound generates other graphs with "vertical" dotted lines be- 
tween the lower and upper half of the circle. To prove our bound inductively 
we need therefore to generalize slightly the class of mirror graphs and their as- 
sociated Gaussian measures to a larger category of graphs G U G U V, called 
generalized mirror graphs or GM graphs and pictured in Figure [51 They are 
identical to mirror graphs except that they can have in addition a certain set V 
of "vertical" dotted lines between the lower and upper half of the circle, again 
without any crossing. 




Figure 6: The generalized mirror graphs 

There is a corresponding measure dvc,v with similar rules; there is a single 
w parameter for each pair of dotted line and its mirror, in particular there is 
a w parameter for each vertical line, Again the covariance between two fields 
belonging to two mirror regions is the infimum of the w parameters crossed from 
one mirror region to the over, staying entirely in e.g. the lower half-disk. The 
upper half-part is still the complex conjugate of the lower half-part. The order 
of a GM graph is again the total number L — 2n + \ V\ of dotted lines and its 
amplitude is given by a pointwise positive integral similar to (|20|) : 
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A GuGuv( x >y) = ^ L / dvGuv j [ Y[ d i x l d i x l ] [ dye] 
J J eec lev 

2n+\V\+l 

Cj((Ti, Zi-l, Zi)Cj((Ti, Zi-l, Zi), (21) 

i=l 

where the z's and z's are either xis, xg's or yes according to the graph. 

Defining the integrand I G uguv i x i v) ofaGM graph so that A GuGuV (x, y) = 
J dv GuV I GuGUV {x,y), we have: 

Lemma 4.2. For any GM graph we have, uniformly in a, x and y: 

I Gu Guv( x ,y)<(K\X\) L M^. (22) 

Inded the quantity Iguguv( x jV) ^ s exactly the same than a pressure graph 
but with two fixed points and some propagators replaced by complex conjugates, 
hence the proof through the norm estimates of Lemma 12.11 is almost identical 
to the one of Theorem 12.21 

We now write the resolvent step which results in an integral Fredholm in- 
equality for the supremum of the amplitudes of any generalized mirror graph. 

Let us define the quantity 

r L (x,y)= sup W L/2 A G[jG[jV {x,y). (23) 

GM graphs G,V | L(G)=L 

We shall prove by induction on L: 

Lemma 4.3. There exists some constant K such that for A small enough 

T L (x, y) < KM^(e- cM ^ x -^ + \X\^ 2 J dze- cMi ^T L (z,y)y (24) 



From that lemma indeed obviously follows 
Lemma 4.4. There exists some constant K such that for A small enough 

F L (x,y) < KM Aj e- cMi \ x - y \. (25) 

Indeed iterating the integral Fredholm equation ([2"3]) leads obviously to (|25l) . 

Taking (|2"Tj) and (|23|) into account to reinstall the \ L I 2 factor, considering 
the equation L — 2n + V and taking a square root because of (p~9|) . Lemma [4~T1 
is then nothing but Lemma I4T41 for the particular case V = 0. 

The rest of this section is therefore devoted to the proof of Lemma POl by a 
simple induction on L. 
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If L = 0, To{x,y) = J duCj(a,x,y,)Cj(a,x,y,). Expanding the Cj(a,x,y) 
propagator, we get 

Fo(x,y)= / du[Cj(x,y) - iV\ dzCj(x,z)a(z)Cj(a,z,y)]Cj(a,x,y). (26) 



For the first term | J dvCj (x, y)Cj (a,x,y)\, we simply use bounds (JTJ) and (|22l) in 
the case L = 0. For the second term we Wick contract the a field (i.e. integrate 
by parts over a). There are two subcases: the Wick contraction hits either 
Cj(a, z, y) or Cj(a, x, y). We then apply the inequality 

\ABC\ < ^(M 2j \B\ 2 + M~ 2j \C\ 2 ), (27) 

which is valid for any positive A. In the first subcase wc take A = J dzCj(x, z), 
B = Cj{a, z, y) and C = Cj(a, z, z)Cj(a, x, y), hence write 



| / dzCj(x, z)Cj(a, z,z)Cj(a, z,y)C J (cr, x,y)\ < 

^CA^z) [ M 2j\ Cj{(7) Zj y) |2 + M~ 2j \Cj(a, z, z)C 3 (a, x, y)\ 2 ] (28) 

and in the second subcase we write similarly 

| J dzCj(x,z)Cj(a,z,y)C j (a,x,z)C j (a,z,y)\ < 

Az ^A [M^\C,(a, z, y)\ 2 + M-^\Cj(a, x, z)Cj(a, z, y)\ 2 ] . (29) 

Using the uniform bound (|22|) on the "trapped loop" \Cj (a, z, z)\ 2 or Cj (a, x, z)\ 2 
in the C term we obtain 

T (x,y) < KM 4j e- cMi \ x -^ + \X\K(r (x,y) 

+M 4 > J dze- cM3 \*- z \T (z,y^j (30) 

so that hence Lemmas and B~4l hold for L = 0. 

We now assume that (|2~4"|) , hence also (|23|) , is true up to order L and we want 
to prove (|24[) at order L + 1. Consider a GM graph of order L + l. If V > 1 we 
can decompose it as a convolution of smaller GM graphs: 



A GuGuv( x >y) = x J d yi A G 1 uG 1 ( x ,yi) A G 2 i>G 2 uv 2 (yi>y) ( 31 ) 

with total orders L\ for G\ and L 2 for G2, V2 = V — {1} strictly smaller than 
L + 1. Applying the induction hypothesis (|25|) to these smaller GM graphs we 
get directly that 

sup \M- {L+1)/2A GuGuv(^y) <KM^e- cM3 \*-y\. (32) 

G,V\L(GUGUV)=L+1,V>0 
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Hence we have now only to prove (M} for mirror graphs with V = 0. Con- 
sider now such a mirror gr aph G. Because of the |A|~ L / 2 in (23]), we should 
remember that we have only a remaining factor |A| L / 2 to use for our bounds on 

Starting at x we simply expand the first resolvent propagator Cj(a, x, x\) as 
Cj(x, xi) — J dzCj(x, z)i^/\a{z)Cj{a, z, x\). 

For the first term we call Xi 1 the point to which x\ is linked by a dotted line 
and apply a Schwarz inequality of the (|27j) type, with: 



A = J dxiCj{x,x\), (33) 
B = dxi Cj(a,x i - 1 ,x i ), 

ii + l<i<2n ii + l<i<2n+l 

/In 
dxi Cj(<r, Xi_i,Xi) JJdxj Cj(cr,Xi-i,Xi). 

2<£<ii — 1 2<i<ii i=l l<i<2ra+l 

It leads, using again the norm bounds of type (2"2")l on the "trapped loop" in the 
first part of C, to a bound 



\\\ 1 / 2 K(r L (x 1 y) + M 4 > J d Xl e- cMi \ x - x -\T r {x u y) 



(34) 



for some r < L. Applying the induction hypothesis concludes to the bound 

(23). 

Finally for the second term we Wick contract again the a field. There are 
again two subcases: the Wick contraction hits either a Cj or a Cj. Let us 
call i the number of half-lines, either on the upper or on the lower circles, which 
are inside the Wick contraction, and x^, ... Xi k or x^, ... Xi k the positions of 
the dotted lines crossed by the Wick contraction. 

We have now two additional difficulties compared to the L = case: 

• we have to sum over where the Wick contraction hits, hence sum over i 
(because the Wick contraction creates a loop, hence potentially dangerous 
combinatoric) . The solution is that the norm bound on the "trapped loop" 
in the C term of (j27|) erases more and more coupling constants as the loop 
gets longer: this easily pays for choosing the Wick contraction. 

• the dotted lines crossed by the Wick contraction should be kept in the A 
term in inequality (27}. In other words they become vertical lines at the 
next step, even if no vertical line was present in the initial graph. This is 
why we had to extend our induction to the category of GM graphs. This 
extension is what solves this difficulty. 

We decompose the amplitude of the graph in the first subcase of Figure [7| as 

/ dzdx il ,...dx ik Cj(x,z)TL Xii! ,„ Xik (z,z)R Xiit ... Xik (z,y)S(x,y) (35) 
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z 



Figure 7: The Wick contraction 

with hopefully straightforward notations, and we apply the Schwarz inequality 
(|27|) . with: 

A = \X\ l/8 J2j dzdx n ,...dx lk JCj(x,z), 

B = Rx il ,...x ik (z,y), 

C = \\r i / s TL Xn ,... Xik (z,z)S(x,y). (36) 

Now the first remark is that z|A| i//8 is bounded by K for small A so we need 
only to find a uniform bound at fixed i. 

The ^4|£?| 2 is a convolution of an explicit propagator bounded by |T]) with a 
new GM graph (with vertical lines which are the crossed lines at Xi t , ...Xj fc ) either 
identical to G or shorter. If it is shorter we apply the induction hypothesis. If 
it is not shorter we obtain a convolution equation term like in the right hand 
side of (gH). 

The A|C| 2 contains a trapped loop TL with i vertices. Each half-vertex 
of the trapped loop has only IAI 1 / 8 because of the |A| -1 / 8 factor in (|36f . The 
trapped loop is again of the GM nature with vertical lines which are the crossed 
lines at x^, ...Xi k . But we can still apply the bound (|22| to this trapped loop. 
Therefore the bound on the sum of the A|i3| 2 and ^4|C| 2 is again of the type 

(IMIl- 
Finally the second subcase, where the Wick contraction 4- hits a Cj, is 
exactly similar, except that the "almost trapped loop" is now something of the 
type TL(x,z) rather than TL(z,z). But the bound (|22|) also covers this case, 
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so that everything goes through. 

Collecting the bounds (|34|) in every case completes the proof of Lemmas 14.31 
and 14.41 for Tl+i. This concludes the proof of Lemmas 14.31 and 14.41 for all L. 

5 Further topics 
5.1 Higher functions 

The analysis of the 2p point functions is similar to that of the previous section. 
The general 2p point function S c (xi, X2 P ) defined by (fT4|) contains p resolvents 
of the Cj(a) type and a certain number of loop vertices joining or decorating 
them. Turning around the tree we can still identify the drawing as a set of 
decorated resolvents joined by local vertices or dotted lines as in Figures |S] and 
O which are the analogs of Figures [3] and [4] This is because any chain of loop 
vertices joining resolvents can be "absorbed" into decorations of one of these 
resolvents. 



xi x 2 




Figure 8: A connected 4 point function 




Figure 9: The "half-disk" representation of that connected 4 point function 

The factor 2pl in (fT5j) can be understood as a first factor 2pU to choose the 
pairing of the points in p resolvents and an other p\ for the choice of the tree 
of connecting loop vertices between them. We can again bound each term of 
the initial expansion by a "mirror" term pointwisc positive in a with p disks as 
shown in Figure [TU1 

A Lemma similar to Lemma 14.11 is again proved by a bound on generalized 
mirror graphs such as those of Figure [10] but with additional vertical lines inside 
the p disks. This bound is proved inductively by a single resolvent step followed 
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Figure 10: The mirror representation of the same connected 4 point function 



by a Fredholm bound similar to Lemmas 14.31 and 14.41 Verifications are left to 
the reader. 

5.2 Large/small Field Expansion 

To prove the tree decay of the 2p-point connected functions as external argu- 
ments are pulled apart, it is possible to replace the Fredholm inequality of the 
previous section by a so-called large/small field expansion. It still relics on a 
resolvent expansion, but integration by parts is replaced by a probabilistic anal- 
ysis over a. We recall only the main idea, as this expansion is explained in detail 
in p~5l[20] but also in a very large number of other earlier publications. 
A lattice V of cubes of side M~i is introduced and the expansion is 

1 = 11 [x(J A M 4j \X\ e a 2 (x)dx) + [1- X ( J A M^\X\^ 2 (x)dx)]j (37) 

where x is a function with compact support independent of j and A. 

The small field region S is the union of all the cubes for which the \ factor has 
been chosen. The complement, called the large field region L, is decomposed 
as the union of connected pieces Lk- Each such connected large field region 
has a small probabilistic factor for each of its cube using e.g. some standard 
Tchebycheff inequality. 

The field is decomposed according to its localization as a = as + J2k a L k - 
Then the resolvent Cj(a,x,y) is simply bounded in norm if x and y belong to 
the same Lk region because the decay is provided by the probabilistic factor 
associated to Lk- 

The as piece is expanded according to resolvent formulas such as 

Cj(as,x,y) = Cj(x,y)-iV\ dzCj(x, z)a s (z)Cj{as, z,y), (38) 
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which can be iterated to infinity because the as field is not integrated with the 
Gaussian measure but bounded with the help of the small field conditions. 

Then inside each connected large field region Lj~ the resolvent Cj{aL k ,x,y) 
is simply bounded in norm. The decay is provided by the probabilistic factor 
associated to L^. Between different connected large field regions, the decay is 
provided by the small field resolvent expansion. 

However one advantage of the loop expansion presented in this paper is to 
avoid the need of any lattice of cubes for cluster/Mayer expansions. If possible, it 
seems better to us to avoid reintroducing a lattice of cubes in such a small/large 
field analysis. 

5.3 Multiscale Expansions 

The result presented in this paper for a single scale model should be extended 
to a multiscale analysis. This means that every loop-vertex or resolvent should 
carry a scale index j which represents the lowest scale which appears in that 
loop or resolvent. Then we know that the forest formula used in this paper 
should be replaced by a so-called "jungle" formula [6] in which links are built 
preferentially between loop vertices and resolvents of highest possible index. 

This jungle formula has to be completed by a "vertical expansion" which 
tests whether connected contributions of higher scales have less or more than 
four external legs of lower scales, see e.g. |15j . A renormalization expansion then 
extracts the local parts of the corresponding two and four point contributions 
and resums them into effective couplings. In this way it should be possible to 
finally complete the program |15j of a Bosonic renormalization-group-rcsummcd 
expansion whose pieces are defined through totally explicit formulas without 
using any induction. Indeed the missing ingredient in |15j . namely an explicit 
formula to insert Mayer expansions between each cluster expansion, would be 
totally avoided. The new multiscale expansion would indeed not require any 
cluster nor Mayer expansion at any stage. 

The expansion would be completed by auxiliary resolvent expansions, either 
with integration by parts in the manner of section [4] or with a small/large field 
analysis as in subsection 15. 21 above. This is necessary to establish scaled spatial 
decay, which in turn is crucial to prove that the renormalizcd two and four point 
contributions are small. But these new auxiliary expansions shall be used only 
to prove the desired bounds, not to define the expansion itself. 

5.4 Vector Models 

The method presented here is especially suited to the treatment of large N 
vector models. Indeed we can decompose a vector </> 4 interaction with an inter- 
mediate scalar field as in ([2|) but in such a way that the flow of vector indices 
occurs within the loop- vertices. Every loop vertex simply carries therefore a 
global N factor where N is the number of colors. Hence we expect that the 
loop expansion presented here is the right tool to glue different regimes of the 
renormalization group governed respectively e.g. in the ultraviolet regime by 
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a small coupling expansion and in the infrared regime by a "non-perturbative" 
large iV expansion of the vector type. This gluing problem occurs in many 
different physical contexts, from mass generation of the two-dimensional Gross- 
Neveu [3U] or non-linear a-model [5T] to the BCS theory of supraconductivity 
[2"2"] . These gluing problems have been considered until now too complicated in 
practice for a rigorous constructive analysis. 

5.5 Matrix models and <f>^ 

The loop expansion is also suited for the treatment of large N matrix models 
and was in fact found for this reason [T7] . Our first goal is to apply it to the full 
construction of non-commutative (j)^ 4 |18j . either in the so-called matrix base 
HHH] or in direct space [2"5] . 

One needs again to develop for that purpose the multiscalc version of the 
expansion and the resolvent bounds analogs to section[4]or subsection [5T2] above. 
Indeed neither the matrix propagator nor the Mehler x space propagator are 
diagonal in the corresponding rcprcscntations/footnotcTherc is an interesting 
exception: the matrix propagator of </>| 4 becomes diagonal in the matrix base 
at the very special ultraviolet fixed point where fl, the Grosse-Wulkenhaar pa- 
rameter, is 1, Of course the general non-diagonal case has to be treated.. 

Ultimately we hope that better understanding the non-commutative models 
of the matrix or quasi-matrix type should be useful in many areas of physics, 
from physics beyond the standard model [571 HE] to more down to earth 
physics such as quark confinement or the quantum Hall effect [30] ■ 
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